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ABSTRACT 

We calculate the predicted UBVRIJHK absolute magnitudes for models of supernova 
progenitors and apply the result to the case of supernova 2005cs. We agree with 
previous results that the initial mass of the star was of low, around 6 to 8 M Q . 
However such stars are thought to go through second dredge-up to become AGB 
stars. We show that had this occurred to the progenitor of 2005cs it would have been 
observed in JHK pre-explosion images. The progenitor was not detected in these 
bands and therefore we conclude that it was not an AGB star. Furthermore if some 
AGB stars do produce supernovae they will have a clear signature in pre-explosion 
near-infrared images. Electron-capture supernovae are thought to occur in AGB stars, 
hence the implication is that 2005cs was not an electron-capture supernova but was 
the collapse of an iron core. 

Key words: stars: evolution - supernova: general - stars: AGB - supernova: 2005cs 
- infrared: stars 



1 INTRODUCTION 

Core-collapse supernovae (SNe) are the spectacular events 
associated with the death of massive stars. They occur once 
the core is no longer supported by nuclear-fusion reactions 
or electron degeneracy-pressure. The core collapses to form 
a neutron star or, if the core is massive enough, a black 
hole. In the neutron star formation a large neutrino flux 
is produced that transfers a large fraction of energy to the 
stellar envelope that is ejected and produces the observed 
luminous display. 

There are two main evolutionary paths that lead to 
core-collapse. The first is the widely known iron core- 
collapse. This occurs in stars initially more massive than 
about 10 MeQ where nuclear burning progresses all the way 
to production of an iron core. Iron being the most stable 
element does not provide energy by further fusion reactions 
so the core collapses. The second core-collapse path is re- 
stricted to stars below about IOMq. In these stars nuclear 
burning progresses no further than carbon burning. After 
carbon burning neutrino cooling reduces the temperature of 
the core preventing further nuclear reactions. Core collapse 
occurs when the oxygen-neon (ONe) core reaches the Chan- 
drasekhar mass (Mch ~ IAMq) and electron degeneracy- 
pressure can no longer support the core. The central den- 
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1 The value varies with the details of convection in stellar models. 
Here we quote masses using convective overshooting. 



sity increases until a density of around 10 9 ' 6 g/cm 3 is 
reached. Then electron-capture by magnesium- 24 ( 24 Mg) 
and/or neon-20 ( 20 Ne) reduces the electron degeneracy- 
pressure further and the collapse accelerat es (|Mivaii et al.l 
ll98d : lGutirrez. Canal fc Garca-Berroll2005l ). 

Electron-capture SNe are thought to occur predomi- 
nately in massive Asymptotic G iant-Branch (AGB) stars 
jWooslev. Heger fc Weaverll2002l ). AGB evolution occurs af- 
ter helium burning in stars from 0.8 to about 8 Mq. Af- 
ter the main sequence and formation of a helium core a 
star will ascend the red giant branch. When this occurs the 
convective envelope penetrates into the helium core mixing 
products from hydrogen burning to the surface, this process 
is known as dredge-up. It decreases the hydrogen, carbon 
and oxygen abundance while increasing helium and nitro- 
gen abundance. After some time helium will ignite in the 
core and the star will move back to the blue. When helium 
burning ends and a carbon-oxygen (CO) core is formed the 
star moves back to the giant branch and dredge-up occurs 
for the second time. This time dredge-up penetrates deeper 
and mixes up nearly all the helium leaving just a thin layer 
covering the CO core and the hydrogen and helium burn- 
ing shells in close proximity. This arrangement is unstable 
and nuclear burning progresses as a series of pulses. The 
hydrogen burning shell burns until the layer of helium is 
thick enough to ignite. This helium rapidly burns until it 
is exhausted and the cycle restarts. AGB stars are more 
luminous than red giants of the same mass because the hy- 
drogen shell is at a much higher temperature in close prox- 
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imity to the helium burning shell. In the most massive AGB 
stars core carbon-burning occurs around the time of 2nd 
dredge-up and so there is an ONe core (carbon burning ex- 
tinguishes after the CO core is completely converted to an 
ONe core and never reignites). This means that if the core 
grows via thermal pulses it may reach the Mch and pro- 
duce an electron-capture SN. It is a race between mass-loss 
and core growth as to whether a SN will o ccur. The evolu - 
tion of these stars has been well stud ied, e.g.lNomotol (| 1987h : 
Ritossa. Garca-Berro fc IbenI dl999f) andlSiessl d2006f) While 
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Eldridge fc ToutJ (|2004l ) and IPoelarends et al.l (2006) have 
concentrated on the mass-ranges over which these events 
occur and when these SNe might occur. 

It is surprising that two quite different mechanisms for 
producing a neutron star could lead to similar type IIP SN 
events. Type IIP SNe are classified by having hydrogen lines 
in their spectra and a long plateau phase in their lightcurves 
during which the luminosity is roughly constant. To produce 
this type of SN a massive and extended hydrogen envelope 
is required. This envelope which produces the SN display 
does not differ substantially between the two paths of core- 
collapse. Discriminating between these two SN types can 
only be achieved by either observing differences in the nucle- 
osynthesis products from the different collapse mechanisms, 
or by observing the progenitor star. In this letter we demon- 
strate the latter method for SN 2005cs. 

During the SN explosion nucleosynthesis occurs and one 
product is n ickel-56 ( 56 Ni) from explosive bu rning of oxygen 
and silicon (IWooslev. Heger fc Weaverll2002l ). 56 Ni provides 
the late time SN luminosity as it decays to cobalt-56 ( 56 Co) 
and then iron- 56 ( 56 Fe). SN 2005cs produced a lower than 
average a mount of nickel in the explosion, of the or der of 
0.01 M (|Pastorello et al.ll2006l ; iTsvetkov et al.ll2006l l. Pas- 
torello et al. (in prep) who use image subtraction to remove 
the background to accurately estimate the sources late-time 
SN luminosity. They suggest it could be as low as 0.004 
Mq. Ty pical nickel masses a re for example 0.016 Mq for SN 
2003gd (|Hendrv et al.ll2005l ) or 0.075 M for SN 1987A. To 
produce nickel, oxygen and silicon are required. Low nickel 
mass therefore indicates that there was very little of these 
two elements around the collapsing core and such a struc- 
ture occurs within massive-AGB stars. Thus the progenitor 
of 2005cs is a good candidate to be such a star. Further- 
more simulations of collapsing ONe white-d warfs indicate 
that they produce less than 0.001 Mq of 56 Ni (|Dessart et al.l 
120061 ). 

In this letter we first discuss our synthetic photometry 
method to calculate UBVRIJHK magnitudes for our stel- 
lar models. We then compare model predictions to the ob- 
servational progenitor limits for SN 2005cs. It was detected 
in a pre-SN / band image but not in deep JHK observa- 
tions. We show that the latter can provide a strong limit on 
whether 2nd dredge-up had occurred 



2 SYNTHETIC UBVRI & JHK MAGNITUDES. 

Stellar evolution codes only produce a few details of the ob- 
servable characteristics such as the luminosity, radius, sur- 
face temperature and composition where the surface is de- 
fined to be where the average optical depth reaches two 
thirds. Comparing models to observed stars can only be 
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Figure 1. Hertzsprung-Russel (HR) diagrams comparing the 
STARS models (red lines) with the Geneva models (blue lines). 
Prom bottom to top the initial masses are 5, 7, 9, 12, 15 
and 20Mq. Upper panel is plotted with log(L/L,Q) versus 
log(T e ff/K) and the lower panel is plotted by the absolute V 
magnitudc and the V-I colour. 



done accurately from spectroscopic observations of the stars. 
However spectroscopic data are generally not available for 
SN progenitors, apart from the rare case of SN 1987A. Re- 
cent work in searching for information on the progenitors 
of supernovae has focused, successfully, on finding observa- 
tions through numerous broad-band filters (combinations of 
UBVRIJHK). Calculating stellar parameters from broad- 
band photometry obviously provides less accurate results 
than from spectroscopic analyses, but the method is the only 
viable one open to us unless a Milky Way, or possibly Local 
Group SN occurs. Previous attempts to determine stellar 
parameters of SN progenitors assigned bolometric correc- 
tions based on stellar type (effective surface temperature) 
for the proge n itors estimated from broad-band colours e.g. 
I Smartt et al.l |2004). However an alternative method is to 
predict the photometric magnitudes of stars along theoreti- 
cal stellar evolution tracks usin g atmosphere models . 

We used the method of iLeieune fc Schaererl (1200 ll ) 
to perform synthetic photo metry on the Cam bridge 
STARS evolution models (|Eldridge et al.l l2006h (see 
http://www.ast.cam.ac.uk/~stars for more details). This 
involves using the BaSeL model atmosphere library 
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Figure 2. The comparison of observed and predicted I — H 
colours for RSG s. The thick black lines are the observed colours 
from [Elias, Frogel & Humphreys] l ll985t ). The red lines are the 
STARS models, while the blue lines are from the Geneva models. 



l|Westera. Leieune fc Buser! Il999l) to work out the flux in 
differ ent broad-band filters jBessell Il99fj| ; iBessell fc Brett] 
1988) for different surface temperatures, gravities and metal- 
licities. Then interpolating in this grid to obtain colours for 
each model. To calculate the absolute magnitudes we used 
a theoretical spectrum of Vega from the same atmosphere 
library, setting all colours to zero. In this we assumed a ra- 
dius of 3.1 Rq for Vega which is larger than its known radius 
of 2.7 R Q (|Ciardi et alj|200ll ). However using the theoretical 
Vega spectrum it provides the correct flux for Vega in the 
V band. The adopted radius is the primary source of error 
in the absolute magnitudes but does not affect the colours. 

The resultant colours agree with the colours from th e 
commonly used Geneva models (jLeieune fc Schaererl l2001). 
Both sets of models are plotted in theoretical Hertzsprung- 
Russel diagrams in Figure Q] The slight disagreement is due 
to differences in the two evolutionary codes, such as up- 
dated opacity tables, different initial abundances and differ- 
ent mass-loss rates. Also our models extend further in time 
as they are calculated to the beginning of neon burning while 
the Geneva models end after core carbon-burning and do not 
have 2nd dredge-up. This means the track endpoints of our 
models are more luminous and slightly cooler. Our models 
for 8 Mq and below experience 2nd dredge-up. 

It is difficult to estimate the uncertainty in our colours. 
The largest errors are due to the assumption of the radius 
for Vega and also assuming that all filters have an absolute 
magnitude of zero for Vega. We estimate the error in the 
absolute magnitudes due to this are around ±0.3 mags. This 
was estimated by using different methods to calibrate the 
zero point of our system such using the colours of Sirius and 
the Sun, and varying the radius of Vega. 

We have checked the predicted colours again st thos e 
of observed RSGs by lElias. Frogel fc Humphreys] (QUI). 
They supply the average colours expected for different 
spectral types of RSGs. However recently there has been 
a reassessment of the effective temperatures of RSGs by 
iLevesque et alj (120051 ) which we take into account in our 
analysis. We show examples of the comparisons in Fig- 
ure [2] The agreement of the m odels with the data of 
lElias Frogel fc Humphreys! (|l985h is poor with the average 



discrepancy of around 0.5 to 0.7 mags for all the colours. 
We have also calculated the expected colours for RSGs 
from the MARCS atm osphere models which were used by 
ILevesque et alj (|2005h to derive their temperature calibra- 
tion. We find that the colours are close to those from the 
atmosphere models we use with a discrepancy of less than 
0.1 mags. 

It is more difficult to compare the predicted colours 
of the AGB stars with observations. While lists o f AGB 
colours exist dKerschbaum. Lebzelter fc Lazardl200ll ). AGB 
stars come from a wide range of initial masses (0.8 < 
M < 7Mq) and therefore have a wide range of core 
masses and luminosities. However the AGB stars with larger 
cores, those closest to the point of core-collapse, tend to 
be the most luminous one. The predicted colours of our 
A GB stars broadly agree with the observ ed AGB colours 
of Kcr schbaum. Lebzelter fc Lazard l|200ll ) and lFluks et al.l 
(1994). Also we have compared th e colours to the mor e de- 
tailed AGB atmosphere models of iGroenewe gcn (2006|) and 
again there is a broad agreement but the range of possible 
colours is much greater in these samples. We will use these 
AGB colours as well as those from our models to determine 
whether the progenitor of SN 2005cs was an AGB star. 

The disagreement between the observed and theoreti- 
cal colours is in general greatest at the end points of the 
stellar models where the RSGs are most luminous and ex- 
tended. This may be due to the way the stellar radius is cal- 
culated in theoretical stellar models. Stellar evolution codes 
use wavelength averaged opacity and the photosphere is de- 
fined to be where the optical depth becomes two thirds. 
However in RSG atmospheres opacity is wavelength depen- 
dent which results in longer wavelength (infra-red) emis- 
sion coming from a smal ler radius than the shorter wave- 
length (visible) emission (I Young et alj 12000) . If this radius 
is not comparable to RSG radii then a systematic error 
is introduced to the predicted colours. However there is 
little we can do to correct the models but we can cal- 
ibrate our predicted colours by introducing an empirical 
correction. Therefore we produced a second set of tracks 
wher e for RSGs we correct the JHK magnitudes to agree 
with lElias. Frogel fc Humphreys! fl985!). The corrections we 
adopt are: J c = J + 0.5, H c = H + 0.5 and K c = K + 0.7. 



3 IMPLICATIONS FOR THE PROGENITOR 
OF SN 2005CS. 

The progenitor of SN2005cs was detected in deep Hubble 
Space Telescope (HST) I band images of M51 taken in Jan- 
uary 2005, a few months before the SN occurred in June 
2005. In BVR and JHK filters the progenitor was not de- 
tected. Two studies exist that use the same HST image 
to calculate I band magnitudes for the progenitor but use 
different near infra-red data. The BVR upper limits indi- 
cate that the progenitor could not have been a blue super- 
giant and was a red supergiant no hotter than a K5Ia type 
(|Li et al.|[200e3 ; iMaund et alj|2005h . However t he JHK lim- 
its als o constrain how cool the progenitor was. iMaund et al.l 
|2005l ) made use of ground-based Gemini observations taken 
in A pril 2005 to pr oduce limits on the JHK magnitudes 
while lLi et aH j|200rj ) used HST NICMOS observations taken 
in June 1998 to produce deeper JH but shallower AT-band 
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Figure 3. The predicted progenitor magnitudes from the STARS 
code, with the correction applied, compared to the observed J 
band detection and JHK limits. The empty dia monds a r e from 
iMaund et all l|2005t ) and the filled diamonds from lLi et al1l|2006l ). 
The shaded region is defin ed by the maximum and minimum 
colours for AGB stars from Grocncwcgcn (2006). The solid lines 
are for RSG models while the dashed lines are for AGB models. 



more models with higher values. Therefore the conclusion 
that the progenitor was not an AGB star is firm. 

For the STARS and Geneva models the progenitor 
would be around 6 to 8 Mq. This o f co urse agre e s wit h 
the conclusions of IMaund etail ([20051 ) and lLi et all (|2006h . 
But now we can limit the amount of 2nd dredge-up that 
occurred in the progenitor. We can conclude that the pro- 
genitor did not experience the extreme dredge-up that pro- 
duces an AGB star. However it is possible that the progen- 
itor did experience a small amount of dredge-up decreas- 
ing the helium core mass slightly, by a few tenths of a so- 
lar inass ; _riris_Jspred2cted in some of the models shown 
bv lPoelarends et al.l (|2006h and this would slightly decrease 
the luminosity of the progenitor meaning we would under- 
estimate its initial mass. Therefore a more robust estimate 
of the final helium core mass for the progenitor is between 
1.7 to 2.2 M Q . 

The important result is that the 2nd dredge-up did not 
lead to an AGB star. If it did it would be clearly detected in 
every JHK image available while agreeing with the I band 
detection. Therefore this sets an upper limit for the masses 
of AGB stars and also AGB stars as SN progenitors of 6 
to 8 Mq. This overlaps with the upper limit for AGB stars 
to produce a white dwarf of around 6.8 to 8.6 Mq from 
iDobbie et all (|2006l ). 



li mits on the progenitor. Both s tudies used the observations 
of lElias. Frogel fe Humphreys! <|l985h to limit the spectral 
type of the progenitor star. 

The non-detections in the JHK bands can be used to 
place very sensitive limits on the luminosity and tempera- 
ture of the progenitor. If we decrease the surface tempera- 
ture of the progenitor while maintaining the I band magni- 
tude we must increase the progenitor's luminosity and the 
synthetic JHK magnitudes become greater than the upper 
limits. This restricts the type of progenitor severely. 

In Figure [3] we plot the I JHK magnitudes of the 
STARS models with the correction. In calculating the 
magnitudes we have used a distan ce modulus of 29.62 
|Feldmeier. Ciardullo fc Jacobv|[l997l ). Over these lines we 
then add the I ba n d det ection and the JHK limits from 
both IMaund et all (120051 ) and iLi et all (|2006|). corrected 
for d ust extinction using an Ay of 0.34 JPastorello et al. 
|2006|) and the extinction law of lCardelli. Clayton fc Mathis 
(1989). Without the correction it is not possible to fit 
the / detection and still stay below the JHK limits un- 
less w«_adopt_the very upper ends of the error bars of 
the ILi et al.l (|2006l ) as their limits. With the correction to 
the JHK colours applied in Figure [3] it is easier to fit 
the / band detection and the near-infrared upper limits. 
This supports our method of correcting our model syn- 
thetic photometry to match the observed JHK colours of 
lElias. Frogel fc Humphreys! l|l985t ). 

The most important feature to notice is that if the pro- 
genitor had gone through second dredge-up to become an 
AGB star, its colours would be as shown by the dashed 
lines and it would have been clearly detected in the near- 
infrared bands. This is made more certain by the shaded re- 
gion s howing the range of AGB colours from iGroenewegenl 
( 2006). Most of the AGB colours are far above the JHK 
limits but those with the lowest values for the JHK mag- 
nitudes have very low mass-loss rates and there are many 



4 DISCUSSION OF UNCERTAINTIES. 

The uncertainties in our estimation of the progenitor mass 
come from three sources. Firstly errors in the observational 
photometry, secondly errors in the synthetic photometry and 

third errors due to the lack of understanding of RSGs. 

The 7- band photomet ric-measurements of lMaund et alj 

l|2005h and lLi et all (|2006l ) differ by 0.3 mags. This is mainly 
due to the different methods employed to remove the flux 
from the nearby bright objects. However we can see in Figure 
0a difference of 0.3 mags in the I band actually makes little 
difference to the conclusions we draw. As discussed above, 
it is the JHK limits that constrains the progenitor more 
severely, particularly cool AGB stars. 

There are four uncertainties in the synthetic photome- 
try. The first two are the assumptions of the Vega zero-point 
for all colours and our assumed Vega radius as mentioned 
above that introduces an error of around 0.3 dex. Third, the 
use of the filter functions for the UBVRIJHK bands; the 
observed I band magnitude was in fact estimated from a 
HST filter which is quite different from the Johnson filters. 
However we have compared the colours predicted by the 
different JHK filter functions and find that the difference 
is typically less than 0.05 dex. We also assumed that the 
synthetic spectra are accurate, however any error in their 
calculation will result in incorrect magnitudes. For example 
the models do not include mass-loss. Another uncertainty 
arises from the use of a wavelength averaged optical depth 
to calculate the stellar radii as discussed above. Other uncer- 
tainties in the stellar evolution code, such as use of mixing 
length theory to describe convection could affect our model 
parameters. However these are accounted for in our correc- 
tion factors. 

Uncertainty is also present due to a number of AGB 
star and RSG enigmas; dust, oscillations and asymmetry. In 
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this case dust is not an important uncertainty. Any dust ex- 
tinction in the I JHK bands is m uch less than the U-band 
l|Cardelli. Clayton fc Matlrislll989h . Because our conclusions 
only depend on these bands they are not affected by dust 
extinction unless a very large reddening (of order A v ~ 10) 
is invoked. But this is inc onsistent with all estimated of red- 
dening towards this SN dMaund et all 120051 ; iLi et all l200rj ; 
IPastorello et al.ll2006l ; lBaronll2006l )~ 

A larger problem is whether the p rogenitor wa s oscil - 
lating or pulsating before the explosion. iHeger at all (|l997l ) 
performed calculations to follow the oscillations a red su- 
pergiant may undergo before explosion. They show large 
changes in the surface details of the star. However there is 
no observational evidence how a star might oscillate before 
a SN. We have included by hand oscillations in our stel- 
lar m odels of a similar magnitude as those of IHeger at all 
(|l997l l We find that while the I band magnitude can vary 
a great deal the JHK magnitudes vary only slightly. 

The last problem of asymmetry is b est shown by obser- 
vations of Betelgeuse (lYoung et al|[2000l) . Not only does this 
work show that the radius of a RSG is wavelength depen- 
dent but we also see that the near infra-red light comes from 
deeper in the star than the visible light. This visible light, 
originating from nearer the surface may be more affected by 
star spots and this would make the star look more luminous 
in the I band than the JHK bands. 

These many different uncertainties do contribute to the 
error on the estimation of the initial mass but do not af- 
fect our conclusion that the progenitor did not experience 
2nd dredge-up and was not an AGB star, arguing against 
the core-collapse mechanism being electron capture for SN 
2005cs. It is worthwhile to compare the progenitor mass 
for this SN t o the u pper mass limit for an AGB star from 
iDobbie et all ([2006) . It is striking that t he masses match 
at aro und 6 Mq. Therefore studies such as lPoelarends et all 
(2006) are required to more deeply understand this impor- 
tant change over in behaviour from an AGB star to a SN 
progenitor. 

In conclusion the progenitor of SN 2005cs was a low 
mass star of between 6 and 8 Mq and may have experienced 
slight 2nd dredge-up but not a large amount of dredge-up 
that would have produced an AGB star. Furthermore if AGB 
stars are the progenitors of some SNe they will leave a clear 
signature in near-infrared pre-explosion images. An attempt 
to put limits on SN progenitors in the near infra-red in a 
systematic way is already underway by our group with deep 
imaging of about 50 galaxies in JHK with VLT, ISAAC, 
Gemini and UKIRT (|Maund et alll2005l ). 
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